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Abstract. In this letter, we revisit the orbit problem, which was studied 
in [11213] . In [3], Kannan and Lipton proved that this problem is decid- 
able in polynomial time. They conjectured that the following problem 
is also decidable: whether the orbit could attach in a given subspace of 
dimension less than four. They further conjectured that this problem can 
be solved in polynomial time when the dimension of the given subspace 
is one. In this paper, we show that this conjecture is valid by establishing 
the connection between this problem and the celebrated Skolem prob- 
lem. In the second part of this paper, we study the approximate orbit 
problem, and show that this problem is decidable except for one case. 

1 Introduction 

A sequential machine [415] is called linear if its "next state" function and "out- 
put" function are linear transformations. Such machines have many important 
applications, for instance, computer control circuitry, automatic error-correction 
circuits, digital communication systems etc. It is very natural to study the ac- 
cessibility problem for linear sequential machines, i.e., decide whether there is an 
input x such that on x the machine, starting in a given state q, goes to a given 
state q' . In [I], Harrison proved that this problem can be reduced to a standard 
linear algebra problem, called the "orbit problem" . 

Definition 1 (Orbit problem) Given an n x n rational matrix A, a rational 
initial vector x and a rational target vector y, decide whether there is some k G N 
such that A k x = y. 

Shank showed that this problem is decidable for n = 2 [2], while the general 
decidability result was obtained by Kannan and Lipton [3], They related this 
problem to another linear algebra problem, the matrix power problem, which is 
the polynomial analog of the well-known discrete logarithm problem for fields of 
prime characteristic [BJ. 
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Definition 2 (Matrix power problem) Given two nxn rational matrices A 
and B, decide whether there exists some k £ N such that A k = B . 

Kannan and Lipton show that the "matrix power problem" for general n is 
decidable in polynomial time, so is the "orbit problem" . A natural generalization 
of the "matrix power problem" is the ABC problem studied by Cai et. at. [8]. 
They proved that it is decidable in polynomial time whether there exists natural 
numbers k and j such that A k B : ' = C for given rational matrices A,B,C with 
AB = BA. 

Some generalized versions of the orbit problem have been addressed in liter- 
ature. For instance, 

Definition 3 (Generalized orbit problem) Given an n x n rational matrix 
A, an initial rational vector x, a vector space which is given by a basis B = 
{v\, ■ ■ ■ ,Vd}, decide whether there exists some k £ N such that A k x £ B, i.e., 
whether there exist ji for 1 < i < d such that 



d is called the parameter of this problem. 

In [3J, it is conjectured that the following "generalized orbit problem" in the 
rational field is decidable for the case d < 4, and for the case rf = 1, it is 
decidable in polynomial time. 

In this paper, we show that the "generalized orbit problem" is related to 
a long-standing open problem, the Skolem problem, which asks whether it is 
decidable that the set of zeros of a linear recurrence is empty. 

Definition 4 (Skolem problem) A sequence of rational (algebraic) numbers 
{u m } (m — 0, 1, • • • ) is called a linear recursive sequence of order d if for all 
m > d we have a relation 



where a\,--- , a<2 are fixed rational (algebraic) numbers which are independent 
of m. Given ai, . . . , a,d and uq, . . . , Ud—i, the Skolem problem asks whether there 
exists some k £ N such that Uk = 0. 

The Skolem problem has been widely studied and it is only known to be decidable 
for d < 5 [7 . Note that the Skolem problem can be rephrased as follows: Given 
a rational nxn matrix A, an n-dimensional rational vector x, and linearly 
independent vectors Vi,V2,--- ,v n -i, does there exist some k £ N such that 
A k x belongs to the (n — l)-dimension vector space spanned by V\, v%, ■ ■ ■ , v n -i? 
Also, the "orbit problem" can be regarded as a special case for an n-dimensional 
space, i.e., a point y. 

Similarly, we can also study the following "generalized matrix power prob- 
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Definition 5 (Generalized matrix power problem) Given n x n rational 
matrices A and B\ , B2 ■ ■ ■ Bd with A being non-singular, decide whether there 
exists some k G N such that A k G span{Bi 1 < i < d}. 

We will also study the "approximate orbit problem" which is formulated as 
follows. 

Definition 6 (Approximate orbit problem) Given a rational n x n matrix 
A, a rational initial vector x, a rational target vector y and a rational number 
5 > 0, decide whether there exists some k G N such that 

\\A k x-y\\ < 5, 

where \ \ ■ \ \ is some given norm of the linear space. 

Throughout the paper, we use V to denote the linear vector space; denotes 
the standard unit basis of V and M n denotes all linear mappings on V. It is 
convenient to define a linear bijection 

vec : M n ^ V <8> V 

by vec(Eij ) — a <g> ej , where Eij is the square matrix with all zeros except a 1 
at location (i,j). 

A square matrix U is called unitary if U'U — UW — I with / being the 
identity matrix. It is well known that any unitary U has the spectral decompo- 
sition, i.e., U — V^DV for some unitary V and diagonal unitary D, stands 
for the complex conjugate of V. 

S(A,x) := {x,Ax,A 2 } is used to denote the orbit defined 

by a matrix A and vector x. 

An algebraic number is a root of a polynomial in Q[x] , the ring of polynomials 
in the variable x over the rationals. We will associate an algebraic number with an 
irreducible polynomial with rational coefficients and a sufficiently good rational 
approximation, which uniquely identify the particular root of the polynomial. A 
root of unity is any complex number that equals 1 when raised to some integer 
power n. 

Structure. The rest of the paper is organized as follows. Section 2 establishes 
a connection between the "generalized orbit problem" and "generalized matrix 
power problem", and connects the later problem to the Skolem problem by 
using the Jordan normal form, then solve the rational case where m — 1,2, 3. In 
section 3, we study the "approximate orbit problem" . Firstly, we gives a fully 
characterization of the set of the limit points of the given orbit. After that, we 
provide a method to solve the "approximate orbit problem" when e docs not 
equal to the distance between the target vector and the set of the limit points, 
which could be well approximated using epsilon net. We conclude the paper in 
section 4. 
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2 Generalized Orbit Problem 

One interesting observation is that the rational Skolem problem of order d + 1 
is not harder than the "generalized orbit problem" with parameter d. 

Lemma 1 The rational Skolem problem of order d+ 1 is polynomially reducible 
to the generalized orbit problem with parameter d. 

Proof. For any rational linear recurrence {u n } of order in polynomial time, 
one can construct a (d+ 1) X (d+ 1) rational matrix A, and two rational vectors 
x, y such that 

u n = y T A n x. 

Observe that u n = Q iff A n x lies in the orthogonal complement of y, a d- 
dimensional subspace. □ 

On the other hand, the "generalized orbit problem" with parameter d is not 
harder than the algebraic Skolem problem of order d + 1, which is stated as 
follows: 

Theorem 1 The generalized orbit problem with parameter d is polynomially 
reducible to the algebraic Skolem problem of order d + 1 for fixed d. 

In particular the generalized orbit problem with parameter m < 3 is decidable 
by invoking the decidability of Skolem problem for order less than 5. 

To prove Theorem 1, we firstly employ the technique introduced in [3] to 
show that the "generalized orbit problem" can be reduced to the "generalized 
matrix power problem" , which in turn can be reduced to the algebraic Skolem 
problem. Interestingly, we also show that the "generalized orbit problem" with 
parameter 1 is decidable in polynomial time. 

2.1 Reduction to the Generalized Matrix Power Problem 

In [3 , it is first showed that "orbit problem" is polynomial-time reducible to the 
"matrix power problem" , here we find that their technique is powerful enough 
to show the following lemma. To make our paper self complete, we will give the 
proof detail. 

Lemma 2 The "generalized orbit problem" is polynomial-time reducible to the 
"generalized matrix power problem". 

Proof. Let 

v = max{Z | x, A x x, A 2 x, • • • , A l x are linear independent} 

and 

C= [x\Ax\---\A v x\. 

Thus, Cisanx (y + 1) matrix. 

Since matrix multiplication and rank finding can be done in polynomial time, 
v and C can be computed in polynomial time for given A and x. There are now 
two cases to consider. 
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1 . v = n — 1. Then C is non-singular and therefore 

d d 

A k x = JiVi A k C = Y2,ii[vi\Avi\--- \A v Vi] 

i=l i=l 
d 

&A k = Y,li[vi\Av i \---\A»v i ]C-\ 
»=i 

Let Bi := [vi\Avi\- • •\A v Vi]C~ 1 , which can be computed in polynomial 
time. Clearly, we have A k — •jiBi where A is non-singular. 

2. v < n — 1. Let 13 = [C|L] be an n X n non-singular matrix, which can 
be constructed in polynomial time. Observe that D A t x = ej+i for i = 
0, 1 • • • , v. Hence, 

d d 
A k x = Y lt v t <s=> D~ 1 A k x = D- 1 ^ jiVi 

i=l i=l 

d 

{D- 1 AD) k ■ D~ 1 x = D^J^j.Vi. 

i=i 

We replace A by D~ 1 AD, x by D~~ 1 x, and Vi by D~ 1 Vi. 
For any i, A % x is an n-dimensional vector with in the last n — v — 1 
coordinates. Thus, we obtain on constraint that X)I=i l% v i has no nonzero 
entry among the last n — v — 1 coordinates. 

The rest of the proof is as same as in [3] to show that the problem can be 
reduced to a (y 4- l)-dimensional problem by deleting all but the first v + 1 
rows and columns of A and all but the first v + 1 entries in x and t>j. Then 
replace A by the the (v + 1) x (y + 1) matrix containing the first v + 1 
rows and columns of A, a; by the first v + 1 components of cc, and by 
the first v + 1 components of Uj. A is non-singular now. Then try to reduce 
the problem to the matrix power problem again until it is done. When it is 
done, the rest is several linear equations and a "generalized matrix power" 
equation. 

□ 

Note that for "generalized matrix power problem" A, B\, B2, ■ ■ ■ , Bd, we can 
assume that Bi e span{I, A,A 2 ,--- ,A d }, otherwise, we only need to deal with 
the matrix space span{Bx, B2 ■ ■ ■ , Bd}Dspan{I, A, A 2 , ■ ■ ■ ,A d } and choose an- 
other basis. 

On the other hand, a "generalized matrix power problem": A k = £? =1 7^ 
is equivalent to (A (g> I) k vec(I) = Xlf=i 7i ' vec(Bi). Thus, n x n matrix power 
problem can be reduced to generalized orbit problem of n 2 dimension. 

2.2 The Reduction 

In this section, we will give the reduction from the matrix power problem to 
the algebraic Skolem problem. Our proof heavily depends on the following result 
obtained by Cai [5]- 
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Lemma 3 The Jordan Normal form of a rational matrix can be computed ex- 
actly in polynomial time. The transformation matrix and its inverse can be ob- 
tained exactly. Here, the input size of the problem is the total lengths of the 
binary representation of all the input entries, and the complexity is measured 
in terms of binary operations. The basis transformation can be applied to any 
commuting matrix. 

Now we are ready to give the proof of Theorem 1. 

Proof, (of Theorem [T]) We first reduce the "generalized orbit problem" with 
parameter d to "generalized matrix power problem" with parameter d, that is 
A k = J2i=i liBi with Bi are linear independent and A being non-singular. 

By Lemma El this problem becomes J k = $^_i7»Ci with J being Jordan 
normal form, which is 

d d 

i=i i=i 

Let = vec(J k ), then the i-th entry of z^ equals Pi(k)X k for some polyno- 
mial Pj(-) and eigenvalue A s . 

Choose d + 1 entries of the above equation, says 

d 

i=l 

where a nonzero vector t lies in the orthogonal complement space of the span 
subspace of qi. Obviously, {ut} is some algebraic linear recurrence of order d+ 1 
if we choose in the following way: first choose those entries Pi(k)X k with 
polynomial Pj(fc) being constant, then Pi(k)X k with polynomial Pi(k) being 
degree-one, and so on. 

This completes the proof of Theorem [T] □ 

Invoking that the algebraic Skolem problem of order less than 5 is decidable, a 
direct consequence of Theorem 1 is: 

Corollary 1 The "generalized orbit problem" and the "generalized matrix power 
problem" with parameter less than 4 "is decidable. 

2.3 Polynomial time algorithm to solve the one-dimensional case 

In this subsection, we will show that there is a polynomial time algorithm to 
solve the generalized orbit problem of dimension one, i.e., d — 1 case. Before 
presenting our algorithm, we give the following result obtained in [3] on which 
our proof heavily depends. 

Lemma 4 JJ|/ Given algebraic numbers a and (3, there is a polynomial algorithm 
to solve the equation a k = f3. Further, if there is a solution, the algorithm finds 
all solution (in a polynomial sited encoding). 
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We can state our result as follows. 

Theorem 2 It is decidable in polynomial time to solve A k x — 72/ for given 
rational A, x and y. 

Proof. According to Theorem 1 , we only need to show that there is a polynomial 
algorithm to solve the equation A k = jB for given rational matrices A, B with 
A being non-singular. 

Invoking LemmaH] assume A — diag(Jk 1 (Xi), Jk 2 (^2), ■ • •, Jki(^i)) is in Jor- 
dan Normal form with </fc s (A s ) being a k s x /c s -Jordan block of eigenvalue A s 
(1 < s < I), now the equation 

A k = 1 B^ J ks (\ s ) k =jB ks , 

where B ks is the diagonal block. 

- Case A: If there is some k s > 1, then on the diagonal we have A^ = 76 and 
on the off diagonal we have kA^ 1 = 7c. Since A is non-singular, 7, b, c ^ 0, 
therefore, n = c\ 8 /b. We simply check directly. 

- Case B: k s = 1 for any s, then on the diagonal we get A* = 76^. Also, we have 
that 7,6fc 7^ 0. Now, we obtain that (j^) k = 5^ for any s. Using Lemma |U 
we know that there is a polynomial algorithm to solve that equation, what's 
left is to check directly. 

This completes the proof of this Theorem. □ 
3 Approximate Orbit Problem 

In this section, we study another generalization of the orbit problem, the "ap- 
proximate orbit problem" , that whether the given orbit could attach some given 
open ball. We show that this problem is decidable except for one case. 
Assume an orbit in a linear space with norm | • 1 1, 

S(A, x) = {x, Ax, A 2 x, ■ ■ ■ A m x, ■■■}, 

a target vector y and a radius S. The problem asks whether there exists some 
k e N such that A k x E B(y,S) with B(x,8) = {z : \\x - z\\ < 5} being the 
open ball of radius 6 and centered at y. 

We introduce the concept of limit points of S(A,x). 

Definition 7 p is a limit point if for any e > 0, there exist at least two different 
integers k% 7^ k-i 6 N such that, the open ball B(p,e) contains A kl x and A k2 x, 
i.e., 

\\p- A kl x\\,\\p - A k2 x\\ < e. 
We write Sl for the set of limit points. 

Assume that A, x,y and S are all rational. Our analysis precedes by three steps: 
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- Step 1, we give a parameterized characterization of Sl- 

- Step 2, we check whether Sl is empty. It is not hard to sec that if Sl — 0, 
the "approximate orbit problem" is easy to solve. Otherwise, more precise 
analysis is needed: we approximate the distance between y and Sl as well 
as possible, where the distance defined as following 

D(y,S L )= inf D(y,q)= inf ||y-q||. 

qeS L qeS L 

- Step 3: Assuming that S ^ D(y, Sl)- We separate the following two cases 
by approximating D{y, Sl)- 

• 5 > D(y, Sl)- One can conclude that there exists some k such that 

D(A k x, y) < S, 

in fact, there are infinitely many such k. 

• S < D(y,SL)- Namely, any limit point of S does not lie in the open 
ball of radius 6 centered at y. To handle this case, we need analysis the 
distance between the limit point of S and non-limit point of S. Then, 
we bound the possible k such that A k x <G B(y, S). 

3.1 Characterization of Sl 

In this subsection, we give the full parameterized characterization of Sl, the set 
of limit points of the orbit S. 

Before studying the general case, we note the following well known lemma 
describes the closure of the matrix group generated by one diagonal unitary, 
where the closure is defined according to the matrix norm induced by the given 
vector norm 1 1 • 1 1 . 

Lemma 5 U — diag(ui,u 2 , ■ ■ • ,u n ) is a diagonal unitary with ^ and m being 
no root of unity for any i ^ j . Then, the closure of {U,U 2 , - ■ ■ , U m , • • • } is the 
set of all diagonal unitaries. 

Remark 1. Apply this lemma on orbit S = (A, x) with A being such diagonal 
unitary, we know that Sl = {Vx : V ranges over all diagonal unitaries }. 

Remark 2. It is direct to see that any unitary U lies in the set of the limit points 
of the group generated by itself. By employing the spectral decomposition of U, 
we only need to deal with the diagonal case, then divide all the eigenvalues into 
disjoint set, and apply the above lemma. 

In general, we have the following result: 

Lemma 6 For given S(A,x) = {x,Ax,A 2 x,---A m x,---}, we have an algo- 
rithm which classifies Sl into the following cases, 

- Case A: Sl is empty, and computes a integer N and a constant c > such 
that and for m > N, one have ||A m a;|| > cm. 
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- Case B: Sl is not empty, and computes a non-singular matrix Q and a 
diagonal matrix C whose diagonal entries are zero or with absolution 1 and 
v such that Sl equals to the closure of Q{Cv, • • • , C m v, ■ ■ ■}. 

Proof. Suppose the Jordan decomposition of A is given by A = PJP- 1 with P 
being a nonsingular matrix, and J being the Jordan normal form of A: 

J = diag(J kl (X 1 ) 7 J k2 (X 2 ), ■ ■ •, J fe ,(Aj)), 

where Jfc s (A s ) is a k s x k s - Jordan block with corresponding eigenvalue A s (1 < 
s<l). 

Let z = P- X x. We have that S{A, x) = {Pz, PJz, • • • , PJ m z, ■ ■ ■ }. As P is 
non-singular, Sl can be obtained by applying P on the limit points of 

S(J,z) = {z,Jz,--- ,J m z,---}. 

To consider the limit points of S' , we first deal with the following orbit with 
Jt(A) being the t x t- Jordan block of eigenvalue A and r ^ 0. 

5(J t (A),r) ={r,J t (A)rv , J t m (A)r, • • • }. 

There are four cases, 

- Case i). |A| > 1. Then <S(J t (A)r) has no finite limit point, as there exist some 
constant a > 1 and integer N such that for any m > N, 

||J™(A)r|| > a m 

- Case ii). |A| < 1. Then S(J t (X),r) has only one limit points as there exist 
some constant < b < 1 and integer N' such that for any m > N' , 

lim ||J t m (-V)r|l < b m . 

- Case iii). |A| = 1 and t = 1. Then J t m (A)r = A m n. 

- Case iv). |A| = 1 and t > 1. Then if r has at least two nonzero entries or the 
only nonzero entry does not lie in the last entry, the points of the set will 
goes to infinity, there exist some constant c > 1 and integer N such that for 
any m > N, 

||J t m (-V)r|l > cm. 

Otherwise, for any m, the only non-zero entry of r is the first entry, r\. Then 
the first entry of the J t m (A)r is \ m r\. 

Applying the above result, we are able to check whether S{J,z) is empty or not. 
If it is Sl is empty, then there are two cases: 

- Case 1. According to case i), we know that there is some a > 1 and a integer 
Ni such that 

\\J m z\\ > a m 

is true for any m > N\; 
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- Case 2. According to case iv), we know that there is some c > and a integer 
N% such that 

\\J m z\ \ > cm 

is true for any to > N?. 

Therefore, for to > max{iVi, N2}, we have 

||J m ;z|| > cm \\PJ m z\\ > dm => \\A m x\\ > c'm 

for to > max{iVi, N2}, where c' = |p-i|| as l-P -1 !! being the induced norm of 
P- 1 . 

Case B: Sl is not empty, invoking observation ii), hi) and iv), we know that Sl 
equals to the set of limit points of 

{Pz, PCz, ■■■ , PC m z, ■■■} = P{z, Cz,--- , C m z, ■■■}, 

where C is obtained by modifying J according to z as follows 



C = 



if |A S | ± 1, 

J' k (X s ) otherwise, 



for each 1 < s < I, where J' k (X s ) is obtained by replacing any element by zero, 
but the only element lies in the (1,1) position. 

Now, C is the wanted diagonal matrix whose non-zero diagonal entries are 
absolution 1 and Q = P. That is, the set of limit point of <S is the closure of 
Q{v,Cv,--- ,C m v}. 

In order to see that, we only need to notice that for any to > 1, C m v 
is a limit point of {v,Cv, ■■■ ,C m v}, where v = z. Thus, the limit point of 
{Cv, • • • , C m v} is the closure of itself. □ 

Combing Lemma [S] and Lemma [Sl we give the parameter description of Sl 
as follows. Firstly, we check whether Sl = 0- If this is not the case, we deal 
with the orbit {v, Cv, • ■ ■ , C m v, ■ ■ ■}■ We consider the equivalence = over the 
non-zero eigenvalues of C: two eigenvalues Ai = A2 if A1/A2 is a root of unity. = 
gives rise to a partition {Si}. Without loss of generality, we assume that 



C 



D 




where D = (BjDj, with Dj being a diagonal unitary matrix whose eigenvalues 
are in Sj. 

There is an integer Nj such that D^ J equals to / up to some global phase. 
Thus, there is a universal N such that equals to I up to some global phase. 
Now we divide {v, Cv, ■ ■ ■ , C m v, ■ ■ ■} into N orbits, the fc-th orbit is 

S{C N ,C k v) = {C k v, C N+k v, • • • , c N+mk v, ■■■}. 
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If there is no j such that the eigenvalues of Sj are unit root, that is D l j ^ I for 
any I, according to Lemma 5, the set of limit points of S(C N ,C k v) is 

S k L = {@ i e i ^D k r v:a j &Tl}. 

If there is j such that the eigenvalues of Sj are unit root, says Si, according to 
Lemma 5, the set of limit points of S(C N ,C k v) is 

S k L = {D k e,>i e ia W k -v.ajG 11}. 

Thus, we can obtain 

S L = \jQ-S k , 

the union of finite parameter described boundary closed sets, a compact set. 
3.2 Approximate the distance between y and Sl 

In this subsection, we give a simple algorithm to approximate the distance be- 
tween y and Sl- 

For convenience, we assume that no there is no j such that the eigenvalues of 
Sj are unit root. And the idea can be directly used to approximate the distance 
for that case with Si contains unit root eigenvalues. 

Firstly, we provide an algorithm which could well approximate the distance 
between y and S k L = {Q ®£ =1 e la W k ■ v : aj that is 

d k := inf / fc (/ii,/i 2 , ■ • ■ ,A»h)> 

with restrictions that \nj \ = 1 for any j, in other words, dk is the infimum of fk 
over the direct product of h unit circle. The function f k is defined on the direct 
product of h unit circle 

/fc(/"i,M2,--- ,f*h) = \\y-Q®^kD k -v\\. 

The function f k is a continues function (actually it is a Lipschitz function) since 
the norm function || • || is continues, and the domain is a compact set, the 
minimum value is achievable. 

We will use epsilon-net to give a well approximation of d k '- For any given 
e > 0, we can choose a set X of finite number of points in the feasible set such 
that for any point of the domain, there is a point of X t such that the distance 
of this two points is less than e. We choose dk(X e ) as the minimal value of fk 
over the finite points set X e as a approximation of d k , that is, 

d k (X e ) := min fk("i, "2, ■ • ■ , v h ). 

Since the function fk is Lipschitz, we know that dk(X e ) is close to dk- To get 
better approximation of dk, we only need to choose smaller e.. 

For any k, we can have well approximation of dk, thus, d(X e ) — min/c dk(X e ) 
could be made sufficient close to D(y,Si)- More precisely, according to the 
Lipschitz property, for any £ > 0, we can choose e such that for any epsilon net 
X f of the direct product of h unit circle, we have 

d(X e )<D(y,S L )+C 
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3.3 Solve the Approximate Orbit Problem 

In this subsection, we attempt to solve the "approximate orbit problem" . At the 
first step, we check whether Sl is an empty set. 

Theorem 3 If Sl = 0, the "approximate orbit problem" is decidable. 

Proof. Sl is empty, we can find K and c > such that for k > K 

\\D{y,A k x)\\ = \\y-A k x\\ > \\A k y\\ - \\y\\ > ck - \\y\\. 

Thus, we only need to verify whether there is some k < max{ e+ jMI , K} that 
whether \\D{y 1 A k x)\\ < e is true in this case. □ 

Theorem 4 If Sl ^ 0, the "approximate orbit problem" is decidable provided 
thate^D(y,S L )- 

Proof. Firstly, we check which the following cases it is: e > D(y, Sl) or e < 
D(y, Sl) by making good enough approximation of D(y, Sl). 

According to the previous section, we can construct a decreasing series Xj 
satisfying < Xj — D(y, Sl) < 2~- 7 holds for any j. Now we have two new series: 
dj = Xj — e and bj = Xj — e — 2~ 3 for j > 0. It is obviously that 

aj > D(y, S L )-e = dj-e- {dj - D(y, S L )) > dj - e - 2~ J = bj. 

Moreover, aj and bj share the limit D{y, Sl) — e. If e > D(y, Sl), we know 
that there exists j such that aj < 0. If e < D{y,Si), there exists j such that 
bj > 0. The idea to distinguish these two cases is: Check whether aj < 0, if the 
answer is yes, e > D(y,SL)', otherwise, check whether bj > 0, if the answer is 
yes, e < D{y,Si)- When the two answers are both no, Check the same thing 
for j + 1. We can always distinguish the two cases, provided e ^ D(y, Sl), We 
summarize this algorithm below. 

- Case 1: e > D(y, Sl)- There exist infinite many k such that D(A k x,y) < e. 
The reason is that the open ball B(y, e) contains one limit point of the orbit, 
it must contains infinite point of the original orbit. 

- Case 2: e < D(y, Sl)- There is a constant lower bound of D(y, Sl) — e by 
obtain bj — r\ > for some j. 

Note that Sl is a compact set, we know that for any k, we can choose q E Sl 
such that D(A k x, q) = D(A k x, S L ), then 

D(y, A k x) > D(y, q) - D(A k x, q) > D(y, S L ) - D(A k x, S L ). 

If k is some integer that we want, i, e., D(y, A k x) < e, we know that 

D(A k x, S L ) > D(y, S L ) - D(y, A k x) > D(y, S L ) - e > r, > 0. 
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Algorithm 1: Bound e — D(y, Sl) 
input : integer n 

output: nonzero 77. (a number between and e — D(y, Sl)) 
real number b <— 0; 
integer j <— 1; 
while j 7^ do 

Construct satisfying < Xj — D(u,Sl) < 2 _J ; 

real number ctj <— a?j — e; 

real number bj 4— Xj — e 2~ J ; 

if a,j < then 
j<-0; 
?7 <- a,-; 

end 

if bj > then 

j<-0; 

7/ <- ^; 
end 
else 

I + l; 

end 
end 

return s 



According to the proof of Lemma HI we can conclude that there exists a 
nonsingular Q such that A k x = QG k r and q — QCr holds for any q £ Sl, 
where 

G=(™),ond C=(o° 



with I? and being diagonal unitary matrices and i? = ©j . (Xj ) with 



I A , j < 1. Due to Lemma [3] and Lemma HI we can choose D — V . Thus, 



D(A k x,S L ) < D(A k x,q) = \\Q(G k - C)r\\ < s\\{G k - C)r\ \ = s\\R k r'\\, 

where s — \ \Q\ \ > and r' is some vector in the smaller dimensional space. 
Notice that there exists < A < 1 and integer K such that for k > K, we 
have 

\\R k r'\ \ < \ k . 



Therefore, 



j] < D(A k x, S L ) < sX k k < - log A 5. 

s 



The rest is to check whether D(y, A k x) < e holds for k < max{— log A 2 , K }. 
The proof of this theorem is completed. □ 
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4 Conclusion 

In this paper, we revisit the orbit problem. The "generalized orbit problem" 
is studied and the open problem left in the paper [3] is solved. In the second 
part, we study the decidability of so called "approximate orbit problem". We 
first provide a complete characterization of the limit points of the given orbit, 
then we demonstrate a method which gives a well approximation of the distance 
between the target vector and the set of limit points. The assumption that the 
radius does not equal to the distance between the target vector and the set of 
limit points plays a central role in our arguments. If this assumption is valid, 
the "approximate orbit problem" can be solved. Without this assumption, this 
problem seems quite hard as the distance (precise value) is difficult to obtain 
even the given vector norm is I 2 norm. 

By removing the assumption with another one, we have an interesting prob- 
lem: For a given rational orbit S(A, x) with Sl being the set of limit points a 
rational target vector y, whether there is k £ N such that |A fc a; — y\\ < D(Sl,v)- 
That is the radius equals to D(Sl,v), the distance between the target vector 
and the set of limit points. 
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